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Abstract. The recent experimental implementation of condensed matter models in 
optical lattices has motivated research on their nonequilibrium behavior. Predictions 
on the dynamics of superconductors following a sudden quench of the pairing 
interaction have been made based on the effective BCS Hamiltonian; however, their 
experimental verification requires the preparation of a suitable excited state of the 
Hubbard model along a twofold constraint: (i) a sufficiently nonadiabatic ramping 
scheme is essential to excite the nonequilibrium dynamics, and (ii) overheating 
beyond the critical temperature of superconductivity must be avoided. For commonly 
discussed interaction ramps there is no clear separation of the corresponding energy 
scales. Here we show that the matching of both conditions is simplified by the intrinsic 
relaxation behavior of ultracold fermionic systems: For the particular example of a 
linear ramp we examine the transient regime of prethermalization [M. Moeckel and 
S. Kehrein, Phys. Rev. Lett. 100, 175702 (2008)] under the crossover from sudden 
to adiabatic switching using Keldysh perturbation theory. A real-time analysis of the 
momentum distribution exhibits a temporal separation of an early energy relaxation 
and its later thermalization by scattering events. For long but finite ramping times this 
separation can be large. In the prethermalization regime the momentum distribution 
resembles a zero temperature Fermi liquid as the energy inserted by the ramp remains 
located in high energy modes. Thus ultracold fermions prove robust to heating which 
simplifies the observation of nonequilibrium BCS dynamics in optical lattices. 
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1. Introduction 

Stimulated by the development of new experimental techniques in the field of ultracold 
atom gases and pump-probe laser spectroscopy, the simulation and implementation of 
established condensed matter model Hamiltonians in novel experimental setups and the 
study of their nonequilibrium behavior has fiorished recently [Il|2[|3l|4[|5l|6l|71|8l|9l [lOl 
[TT1[12]. In these systems precise and rapid control over parameters can be combined with 
a high temporal resolution for the observation of the many-body quantum dynamics. 
Thus, the switch of an interaction and the analysis of the subsequent response allows to 
address the interplay of particle correlations imposed by interactions and nonequilibrium 
initial conditions. 

For many systems a slow adiabatic change of a two-particle interaction may modify 
the ground state(s) of a system but no excitations to higher energy states occur. 
Therefore, adiabatic switching procedures allow to develop equilibrium properties of 
interacting systems from those of noninteracting ones. Due to the broad applicability of 
this approach they have been counted among the "basic notions" of condensed matter 
theory p^. The most prominent example is Landau's theory of a Fermi liquid [14j: 
Based on its main prerequisite, the adiabatic increase of the fermion interaction starting 
from the noninteracting Fermi gas. Landau could establish a one-to-one relation between 
the eigenstates of the noninteracting Fermi gas and approximate eigenstates of the 
interacting Fermi liquid. This link allows to carry over an approximate particle-like 
description to the interacting system and leads to the celebrated and imaginative concept 
of quasiparticles: reminiscent of the constituent physical fermions, albeit modified 
by a renormalization of their energies and masses, they represent the elementary 
excitations of the interacting system. Note that the quality of this approximation, 
i.e. the considerable time stability of low energy quasiparticles, is owed to phase space 
restrictions for elastic scattering of fermions around the Fermi energy: these restrictions, 
which emerge from the Pauli exclusion principle for fermions, prevent a fast decay of 
approximate quasiparticles into the true eigenstates of the interacting system and are the 
root of the universal applicability of Fermi liquid theory to many interacting fermionic 
systems. 

However, nonadiabatic switching, e.g. of the strength of trapping potentials or 
magnetic fields, of optical lattices or by applying ultrashort laser (pump) pulses may 
drive a system away even from its interacting ground state. Hence, fast parameter 
switching has shown to be a helpful tool to impose nonequilibrium initial conditions onto 
many-particle quantum systems. The limiting case of an instantaneous switch, called 
a quantum quench, provides a convenient starting point for theoretical investigations. 
Afterwards, the evolution of the excited many-body state can be followed in real time 
and diflFerent stages of the relaxation dynamics have been observed. 

In fermionic systems, nonadiabatic switching of the interaction challenges the 
main prerequisite of Fermi liquid theory. Nonetheless, for not too strong sudden 
interaction quenches within the Fermi liquid phase of the Hubbard model, it has 
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been predicted analytically [151 US] , confirmed numerically [H] and furtheron discussed 
pn [T8] that during the initial stage of the relaxation a nonequilibrium description builds 
up which is strongly reminiscent of Fermi liquid theory. The same fermionic phase space 
restrictions which ensure the persistence of quasiparticles in equilibrium restrict a rapid 
relaxation of the momentum distribution in nonequilibrium. Hence, a transient regime 
of prethermalization has been observed. Whereas, for instance, the total kinetic energy 
has already relaxed to its final value, the momentum distribution combines different 
features: a sharp step at the Fermi energy indicates zero temperature while deviations 
from its equilibrium shape represent the nonequilibrium nature of the transient state. 
Hence, during the first phase of the relaxation dynamics no substantial conversion of 
the initial excitation energy into heat can be observed and thermalization is deferred to 
a later time scale [15j. 

In this paper we generalize these observations to finite switching times, assuming 
a linear ramping procedure. This introduces a new time scale, the ramp-up time which 
is defined as the time duration on which the parameter change occurs. This ramp up 
scale can be used deliberately by finite-time 'ramp-up' protocols to prepare a system in 
particular (excited) states which are different from those reachable by sudden switching 
[19]. For instance, ramping on finite time scales has been suggested as a way to access 
the nonequilibrium regime of systems which are sensitive to heating effects since slower 
parameter changes reduce the excitation energy inserted into a system. 

In particular, this argument has been applied for the analysis of superconducting 
BCS systems in nonequilibrium. Theoretical studies have examined an interaction 
quench in the BCS Hamiltonian and predicted interesting nonequilibrium dynamics, 
e.g. oscillations of the order parameter in time [2Ql EJ, 122 However, the 

BCS interaction emerges from a local two-particle interaction between fermions which 
is the tunable switching parameter in experiments with ultracold fermions loaded onto 
an optical lattice. A realistic modelling of nonequilibrium dynamics in ultracold atomic 
gases therefore has to start with a quench of a spatially local interaction. Hence, BCS 
quench dynamics must be discussed against the background of a quenched Fermi liquid. 
Obviously, BCS dynamics can only be expected for temperatures below the critical 
temperature Tc of superconductivity. A first approach based on a golden rule argument 
suggested that overheating can be avoided as energy scales of heating and of the 
nonadiabaticity requirement imposed on the ramp should separate ^6]. However, we will 
show in the sequel that this argument neglects contributions of high energy momentum 
modes which are relevant in commonly discussed ramping models. Therefore, the 
argument turns out inconclusive in the case of linear ramping. Fortunately, the problem 
of overheating is mitigated by two aspects: In another article published in this focus 
issue M. Eckstein and M. KoUar discuss the dependence of the excitation energy inserted 
by interaction ramping on the functional shape of the ramp [27]; this opens space for 
optimizing the ramping procedure. 

We address, on the other hand, the intrinsic robustness of fermionic systems at 
zero temperature to heating. Even for ramping scenarios which do not lead to a clear 
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energy separation of the relevant BCS energy scale and the effective temperature after 
heating we find a temporal separation of the characteristic time scales of BCS dynamics 
and thermalization, respectively. Our observation mirrors the notorious problem of 
cooling a gas of cold fermionic atoms to even lower temperatures and views a current 
technical limitation as a desirable feature for future experiments: On the one hand, 
the Pauli principle suppresses s-wave scattering between the fermions. This hinders 
the equilibration of the gas at a lower temperature after high energy particles have 
been lost by evaporation and poses a great challenge to experimentalists on their way 
to reach a 'zero temperature' Fermi gas. On the other hand, we will show that the 
same mechanism prevents a rapid thermalization of excitations of high energy modes 
in a fermionic system which has been previously cooled down to zero temperature. 
Then a long-lasting transient state of prethermalization emerges during which no 
heating occurs. This behavior of cold fermionic many-particle systems will simplify 
the observation of nonequilibrium BCS dynamics during a long transient time window 
before thermalization sets in; since it roots in Fermi degeneracy caused by the Pauli 
principle it is independent of any details of the applied ramp. 

In this work we perform a careful investigation of the crossover from adiabatic 
to instantaneous linear switching procedures on the grounds of second order Keldysh 
perturbation theory. We follow the evolution of the momentum distribution function 
for original fermions [|] after a linear interaction ramp of finite duration tR and discuss 
the kinetic energy in a long-time limit of the second order result. We observe that the 
corresponding effective temperature is of the same order of magnitude as the critical 
temperature Tc. Hence overheating cannot be excluded by an argument based on energy 
scale separation. However, we still observe prethermalization behavior of the Fermi gas; 
therefore this temperature is only reached after a long transient regime. Since in this 
regime the momentum distribution resembles closely a zero-temperature Fermi liquid we 
conclude that the observation of nonequilibrium BCS dynamics can be expected before 
overheating wipes out all BCS signatures. 

2. Hubbard model with time dependent interaction 

For an analysis of many-body quantum dynamics we implement Fermi liquid theory 
microscopically by the Hubbard Hamiltonian = i?™ + It provides a 

paradigmatic description for correlated fermions on a lattice; it includes spin and 
implements the Pauli principle by restricting to a local state space of dimension 
four. The itinerant properties of the fermions are contained in the kinetic part of the 
Hamiltonian and described as coherent hopping processes between neighboring lattice 
sites with a hopping matrix element h. The chemical potential /x is fixed to ensure half 

:|: Note that we analyze the nonequilibrium properties of the interacting system in terms of the 
unrenormalized eigenmodes of the noninteracting system. This differs from other approaches which 
assume the existence of a quasiparticle picture and discuss directly the behavior of nonequilibrium 
quasiparticle distributions. 
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filling. 

{i,j)a j /cG/C,aG{t,;} 

A momentum space representation is obtained by a Fourier transform. The dispersion 
relation Ck depends on the details of the lattice. For convenience, explicit numerical 
calculations of the momentum distribution are performed on a hypercubic lattice in the 
limit of infinite dimensions {d oc) for which a Gaussian density of states is applicable 
[28]. 

p(e) = exp (-(e//i*)V2) /V2^h' (2) 

/i* is linked to the hopping matrix element by dimensional scaling h — > [28j. For 
convenience, we set h = 1 and fc^ = 1. 

Correlations between the lattice fermions are modeled by a local two-particle 
interaction U acting between fermions of different spin with nja- denoting the local 
occupation 

HT = m (n,t - (n,, - (3) 

Since our interest is in Fermi liquid properties of the model, we restrict to a regime of 
weak interaction strength {U{t) <C Vt). This rules out the observation of any features 
of the Mott-Hubbard phase transition. In consequence, the hopping matrix element is 
a natural choice for the energy unit (/i* = 1). 

To study switch-on processes we assume that the bare interaction U{t) is a time 
dependent parameter. For simplicity, linear ramps of the interaction are considered 
which are characterized by a ramp-up time scale t^. 

{0 t<0 
t/tR 0<t<tR (4) 
1 t>tR 

In the sequel we discuss the effects of variations of the parameter tji onto the real-time 
evolution of the momentum distribution. We assess the dependence of the heating on 
the ramp-up time tR by analyzing the ti^-dependence of the kinetic energy for deviations 
from adiabatic switching procedures. 

3. Keldysh perturbation theory for time dependent Hubbard model 

Perturbation theory for the Hubbard model has been set up previously by various 
authors [29j. Here a real-time Keldysh approach following Rammer allows to 
extract the time dependent behavior of the momentum distribution function N^^^{t) 
from the Keldysh component Greens function {x^t^Gz^x^t^Gz)' Both functions are 
evaluated in second order of the time dependent local two-particle interaction and in 
the initial representation of noninter acting fermions. The later implies that interaction 
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effects become visible only by a redistribution of occupation among the unrenormalized 
momentum modes. 

Nk^^{t) = 5^ / dxe'^'^G^ix, t, a„ x, t, a,) (5) 

Since first-order diagrams vanish due to symmetry, only a single second order diagram 
has to be evaluated for a consistent second order result of the momentum distribution 



function. Details of the calculation can be found in Appendix A or in pT] . 

The second order result describes the dynamics of the momentum distribution 
during a first stage of its relaxation and exhibits the early relaxation of the kinetic 
energy. Higher order corrections, however, induce different behavior on a second time 
scale. We describe this second stage by means of the quantum Boltzmann equation, 
taking the formal long-time limit of the second order Keldysh perturbative calculation 
as an initial condition. It will lead to a thermalization of the momentum distribution. In 
section [6] we will justify this treatment, explain the separation of both relaxation scales 
as a consequence of fermionic phase space restrictions and discuss the characteristics of 
the transient regime of prethermalization. 

4. Evaluation of the momentum distribution 

We have evaluated the time dependent momentum distribution function of original 
fermions in the regime t > tji for arbitrary ramp-up times t^. Within the validity of 
second order time-dependent perturbation theory, i.e. for not too large times t, the 
momentum distribution function reads 

^N^^M) - -U^ j dE ME;n) {^J^^^^^^^^ + 



{E - e,y 

2 sin((^ - ek)tR) cos((E - ek)t) + sin((^ - ek)t){l - cos{{E - 6^)^^)) 



(6) 



tR {E - 

Jk{E]n) represents a fermionic phase space factor which is discussed in more detail 
in [Appendix B[ Here we note that in the case of e^^ > it is nonvanishing only 



for negative energies in the case of < for positive ones. Moreover, it 
asymptotically approaches a quadratic energy dependence for small absolute values of 
i.e. lim^^o J{E]n) ^ E'^. Together, these observations imply that all formal energy 
divergences in equation ^ are cut off or regularized by the phase space factor. 

4.1. No secular terms 

A second glance at the result ^ shows that -contrary to the Keldysh Greens function 



itself (see appendix Appendix A)- the momentum distribution function does not exhibit 
secular terms which are proportional to time t and typically would restrict the validity 
of a perturbative approach to short time regimes. This advantage is a consequence of 
the simplicity of the number operator and of the linear ramping assumption Q. 



Crossover from adiabatic to sudden interaction quenches in the Hubbard model 



7 



4.2. Limiting cases for the linear ramping 

For the correction to the momentum distribution two hmiting cases have been previously 
discussed based on a flow equation calculation [151 ES]- Both limits are reproduced by 
Keldysh perturbation theory. 

4>2.1. Adiabatic limit The adiabatic limit is given by an arbitrary slow linear increase 
of the interaction strength and corresponds to ^ oc. Then only the ti^-independent 
term [E — ek)~^ contributes. 



for a microscopic implementation of Landau's Fermi liquid theory. 

4-2.2. Quench limit The limit of a quantum quench, i.e. a sudden switch-on of the 
interaction, is obtained for tji ^ and represents the maximally nonequilibrium case 
for any linear ramping procedure. Replacing 1 — cos((£' — ek)tR) = 2sin^((£' — ek)T/2) 
the correction to the momentum distribution shows a nontrivial time evolution which 
incorporates the initial build-up of correlations. 



The nonequilibrium nature of this distribution is best analyzed in the long-time limit 
of this second order result. It is deflned for any time dependent function f{t) by 



Note that the notation f{t) intends to resemble the generic time dependence of the 
function f{t) but not a possible time dependence of the long-time average. In this limit a 
characteristic mismatch fi{ek) between nonequilibrium distribution and a corresponding 
equilibrium one -computed for the same interaction strength U- has been found and 
will be discussed in great detail later. 

5. Discussion of the crossover from instantaneous to adiabatic switching 

We discuss the crossover from instantaneous to adiabatic switching by analyzing two 
observables which show signiflcantly different behavior, namely 

1. the time dependent momentum distribution function N^^^{t^tR) and 

2. the kinetic energy E'kin in a formal long-time limit of a perturbative result 

The observed differences illustrate and resolve prior confusion on the relevant energy 
scales of fermionic systems in nonequilibrium. 



hm AiVr^(t,t^) = -U' [dE UE-n) ^ =: AiVf^^ (7) 



This is a truly stationary state. It equals the second order equilibrium result AN^^^ 
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5.1. Time dependent momentum distribution function 

After the interaction has been switched on, correlations between particles develop. In 
the following we observe their impact on the momentum distribution, leading to a 
redistribution of occupation among the momentum modes with time. Features of the 
long-time limit of this dynamics are discussed in a second step. 



5.1.1. Correlation buildup for arbitrary ramping times. The initial buildup of 
correlations is mirrored by a reshape of the momentum distribution. In figure [T] we 
picture this process for five time steps, comparing three ramping scenarios: (1) a fast 
ramp-up of the interaction (black line for ti? = 0), (2) an intermediate ramping speed 
(red line, tji = 4) and (3) the adiabatic limit {tji oc). In the first case the evolution 
of the momentum distribution can be clearly followed: For small times, 1/t-oscillations 
depict a rapid redistribution of the initial occupation (at t = 1) into many high-energy 
modes. With increasing time, occupation is transferred back towards the Fermi energy 
such that in the long-time limit changes to the momentum distribution are most relevant 
around the Fermi energy. 

For intermediate switching times, the process of redistribution of momentum mode 
occupation occurs in a decelerated way. Since we have assumed t > tji in the derivation 
of ([g]) the second time slice {tR = 4, t = 2) does not include possible transient corrections 
which might change the short- time dynamics. 

A comparison of all three ramp-up scenarios in the long-time limit shows that the 
absolute change in the momentum distribution is reduced with increasing ramp-up time 
tR. Most clearly, this is seen at the Fermi energy. 



5.1.2. Mismatch function. A quantitative measure for the deviations of the long time 
limit of the nonequilibrium momentum distribution (tji < oc) from the equilibrium one 
{tji ^ oc) is given by the mismatch function 



While for adiabatic switching on of the interaction lii{ek^tR oc) = 1 is trivial for 
all values of e/^, for finite ramp-up times deviations become visible around the Fermi 
energy. Within a finite environment of approximate radius cr ^ ^/tR the values of 
l^i^k-itR) increases up to a maximum value of two which is reached in the limiting case 
of a sudden interaction quench. This confirms previous observations based on a fiow 
equations approach [16] . 

Fig. [2] shows numerical computations of //(e/^) for different ramp-up times. While 
for strictly sudden switching the maximum value of two is assumed for all energies, 
even for very fast switching {tR = 0.1) a reduction of this value at the band edges 
(bandwidth D=6) is noticeable. On the other hand, even slight deviations from adiabatic 
switching leave an observable 'spike' around the Fermi energy in the mismatch function 
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Figure 1. Build-up of an interacting Fermi liquid description after fast {Ir = 1), 
intermediate {tR = 4) and adiabatic {tR oo) ramping. 



l^{^k)tR ^ !)• This proofs that around the Fermi edge the momentum distribution is 
highly sensitive even for small departures from the adiabatic limit. 

Deviations from the adiabatic behavior are most pronounced at the Fermi energy. 
Hence the (numerical resolvable) peak height of the mismatch function at the Fermi 
energy /x(6i?,ti^) defines a simple measure for the adiabatic character of the momentum 
distribution. Fig [3] displays its dependence on the ramp-up time for some values of tR. 

5.1.3. Intrinsic energy scales of the momentum distribution. From the peak width 
in figure [2] one may be tempted to read off the appearance of a new energy scale 
cr ^ 1/tR m the nonequilibrium momentum distribution. It describes the extension 
of an environment around the Fermi energy 6i? = where significant deviations from 
equilibrium can be observed. Note that it is not equivalent to a scale set by an effective 
temperature since the momentum distribution still exhibits a discontinuity at the Fermi 
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Figure 2. The mismatch function /i(e) is plotted for different ramp- up times. 
Deviations from the value /i(e) = 1 depict the nonequilibrium nature of the transient 
state which is described by the long-time limit of a second order perturbative 
calculation. Those are bounded by an upper value of ii^ep) = 2 in the limit of a 
quantum quench. For small corrections to the adiabatic limit, i.e. for long ramp-up 
times, they are sharply peaked around the Fermi energy. 



energy which indicates zero-temperature behavior. 

The same scale can be extracted from a formal analysis of the momentum 
distribution. Since the sinusoidal functions in ^ can be discussed in two limiting 
cases known as the 'rapid oscillation' and the 'small angle' regime, two energy scales 
= l/tji and 1/t emerge. However, only the first one is relevant in the long-time limit. 
This motivates the discussion of two extremal cases of the internal energy integration 
in ([g]): \E — Ckl ^ l/tR (di) and \E — ek\ <C 1/tR (b). Due to the properties of the 
phase space factor Jk{E^ n) (cf. [Appendix B ) these restrictions give rise to an analogous 



separation of regimes for the external energies ej. ^ l/^i?. 

In the first case, the rapid oscillations and the small 2/t\{E — contributions to 
the integral can be neglected. 



\E-eu\'>l/tR 



= -u' 



I 



dE 



Jk{E,n) 



^^EQU 



\E-ek\->l/tR 



(9) 
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Figure 3. Numerically calculated values of the mismatch function at the Fermi 
energy (circles) are given for different ramp- up times Ir. Within a wide range 
of ramping times, correction to the adiabatic limit {tR oo) can be fitted by 
l^{^F,tR) ^ 1.03 + l.761/tR - 0.921/^2^ + Omi/t% (x^ ^ 10-^, 9 data points were 
fitted); this agrees with the asymptotic expectation /J.{eF^tR oo) = 1. Close to 
a sudden switch there is saturation at the value /J.{eF^tR ^ 0) = 2. The crossover 
between both regimes can be estimated as ^ l//i* = 1 where the inverse power-law 
fitting has broken down. 



For finite ramp-up times and large values of \ek\ the condition (a) is satisfied for all 
values of E. Hence far away from the Fermi energy the nonequilibrium momentum 
distribution equals the adiabatic one ([7|. 

In the second case we consider the long-time limit such that the time dependent 
arguments average out. Since now \E — ek\ < 1/tR the main contributions to the energy 
integral in ([g]) originate from the term 



dEJk{E,n) 



2 (l-cos((£;-e,)t^)) 



/2 



+ 



\E-ek\<l/tR 



Jk{E,n) 



In the case of rapid switching the condition \E — ek\ < l/tR can be satisfied both for a 
large range of energies Ck and, at the Fermi energy, for a large energy interval in E. This 
implies substantial deviations from the equilibrium behavior in an environment around 
the Fermi energy. There it leads to the increase of the mismatch function up to a value 
of two. 

For slower ramping times the restriction cuts off the internal energy integration such 
that the maximum value of two is not reached any more. 
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6. Effective heating for deviations from adiabatic switching 

We start the discussion of the kinetic energy by considering smaU nonadiabatic 
corrections to an adiabatic switching procedure, i.e. we focus on large ramp-up times. 
While the energy scale = l/tji allows to discuss energetically well-separated limiting 
cases of the momentum distribution, it leads to wrong expectations when it is carried 
over naively to the kinetic energy. We will discuss this observation and its consequences 
in the following sections. 

For the purpose of ramp-up dynamics, the kinetic energy AE^^^{tR) is defined as 
the excess of kinetic energy with respect to the corresponding interacting ground state. 



dck Ck 



(10) 



Note that this quantity provides an integrative view on the implications of non-adiabatic 
switching; it includes energy contributions both from the pronounced nonequilibrium 
regime \ek\ <C 1/tR and from a large crossover region \ek\ ^ l/tR of higher energy modes. 



6.1. Full numerical solution of the kinetic energy 



A numerical integration of (10) after insertion of ([g]) exhibits the dependence of the 
kinetic energy AE^^^ (tR) on the ramp-up time (see figure |4]). For small corrections to 
the adiabatic limit, i.e. for long ramp-up times, this behavior agrees with expectations 
formulated on thermodynamic grounds [32j. 



IP 



(11) 



The quadratic behavior can be attributed to the contributions of high energy modes 
with energies larger than a cut-off set by l/tR. This can be seen if one deliberately cuts 
off these contributions in the numerical result: then a different dependence of the kinetic 
energy on the ramping AE^^^ (tR) ^ U^^r^ is rediscovered which has been published 
previously [26l[27j. This indicates that the excitation energy per momentum mode CkSuk 
remains of order of the average AE™^ /D even for \ek\ ^ l/^i? where D represents the 
finite bandwidth. The additional power 1/tR in AE™^{tR) then merely refiects the 
arbitrarily chosen integration limits introduced by the cut-off. 

In the opposite limit of a rapid switch the kinetic energy saturates at the value 
found for the quantum quench scenario (cf. figure [5|). 



6.2. Permits golden rule analysis for a ramping model 

The same observation can be made analytically in the ramping model defined by 
^{t-itR) = AoB(0)(l — e~^/^^) which has been discussed before [26j. In the following 
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Q — o Numerical data 
--■ E(tJ = 0.03t / 



E(tj^) = 0.34tj^ 




Figure 4. The kinetic energy is evaluated numerically for large ramping times (circles). 
Fits for different power-law behavior are given. Obviously, only the inverse quadratic 
fit AE^^^{tR) - agrees weh with the data. 



we will calculate the excitation energy according to Fermi's golden rule. 



0<u;<A 



^nin2(l - n3)(l - n4)5i234 



1234 



(12) 



N{uj) 



We denote the phase space factor of two-particle scattering in the presence of a filled 
Fermi sea by N{(jj) and analyze the role of an energy cut-off A. 

Fourier transforming the time dependent interaction reads X{(jj^tR) = iXo/(jj{l — 
iuotpi) or \\{uj^tR)\^ = |Ao|^/cc;^(l + cc;^t^). In a small frequency limit, the phase space 
factor for a Fermi liquid can be approximated by N{(jj) ^ {Pf/p%)^'^ such that 

A^-^^ r^c.^,^4c^^^4¥ r^M.) ^^^^^^ 



Pf I , ,2 AtR - arctan(AtR) 

7, ^ 
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Figure 5. The kinetic energy in the hmit of fast ramp- up times saturates for a sudden 
quantum quench. 

Now one easily observes that the dependence of the excitation energy on the ramping 
time is a direct consequence of the apphed cut-off. Firstly, the pronounced low energy 
approach taken by Barankov and Levitov cuts off all contributions on the scale A = l/t^^, 
i.e. assuming that for all relevant energy modes utR <C 1. This corresponds to 
an approximation of the Fourier transform by its low frequency asymptotic behavior 
X{(jU^tR) ^ lim^^o A(cl;, ti^) = iX^/uj. Since this approximation mirrors the Fourier 
transform of a sudden switch J^[Q{t)]{uj) = i/u) all details of the ramping model are 
lost. The resulting power-law behavior is AE^^^{tR) ^ . 

If one, on the other hand, one chooses a cut-off which is independent of the ramping 
time (A ^ A(l/ti^)) but small enough such that the low frequency approximation of 
the phase space still holds (this is possible for small corrections to adiabatic ramping) 
one observes the different behavior /\E^^^{tR) ^ . This corresponds to replacing 
the Fourier transform of the switching model by its high frequency asymptotic limit 
lim^^oo A(cl;, ti^) = —Xo/uj'^tR and includes excitations which are characteristic for 
the particular ramping scheme. This approximation agrees well with our numerical 
calculations for linear ramping and is in full agreement with the general analysis provided 
in [27] . Since finite time ramping remains in all cases linked to a continuous increase 
of the interaction strength on the chosen time scale we expect that high frequency 
excitations will be present in many experimental implementations. However, in the 
following we will show that the influence of high frequency excitations on the low energy 
behavior can be delayed by a relaxation bottleneck. Independent of the applied ramping 
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procedure the latter accounts for the emergence of a transient regime characterized by 
prethermahzation. This defers the problem of overheating to late times of the relaxation 
dynamics. 



6.3. Prethermahzation for arbitrary switching times 

As we have already mentioned, the observation of a prethermahzation regime has been 
first reported in the case of a sudden interaction quench in the Hubbard model. Here we 
have discussed the behavior of the momentum distribution function and of the kinetic 
energy on a first time scale set by second order perturbation theory tpx ^ PfU~^ . For 
this time regime and in the case of finite switching times tji <C tpx we have observed that 
the additional kinetic energy remains smaller than the one inserted under a quantum 
quench; moreover, that the momentum distribution function retains a sharp step and 
differs only quantitatively but not qualitatively from that one of the quantum quench. 
In consequence, we note that the phase space restrictions for scattering processes are not 
overcome by finite ramping times. On the contrary, the nonequilibrium character of the 
momentum distribution is less pronounced which makes those restrictions more severe. 
This can be seen when the scattering integral Xk[N^^ {t)] of a quantum Boltzmann 
equation 

^-^=UNk{t)] (14) 

is discussed. Note that the scattering integral vanishes for Fermi-Dirac distribution 
functions at arbitrary temperature and that the quantum Boltzmann equation provides 
an effective description for the evolution of a quasiparticle distribution function. 
Quasiparticles are approximately noninteracting degrees of freedom which absorb most 
of the interaction effects present among the physical fermions. Therefore we construct 
a rough approximation to the quasiparticle mapping by the linear extrapolation 

which strictly holds in the adiabatic limit. In the nonequilibrium case it provides an 
approximate representation of the quasiparticle momentum distribution. 

^^NEQ:QP ^ /^^NEQ _ ^^EQU 



Note that AN^^^'^^ - 0{U^). Since around the Fermi energy |AA^f^^| > |AA^f^^| 
the second order perturbative long-time limit of the quasiparticle distribution shows a 
reduction of the quasiparticle residue and open phase space for quasiparticle scattering. 
Hence the scattering integral, initialized with this long-time limiting state, does not 
vanish but gives rise to a subsequent Boltzmann dynamics of the momentum distribution 



(14). Since its fixed point is a thermal state, thermalization can be predicted. The 
corresponding time scale can be extracted from a linearization of the scattering integral 
w.r.t. AA^^^^'^^ around the zero-temperature initial state. At the Fermi energy the 
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linearized scattering integral can be expressed in terms of the long-time limit of the 
kinetic energy AE^^^ (10) after the initial energy relaxation has been accomplished. [§] 



Q roo o 



Inserting this into (14) and using ^for not too strong derivations from the adiabatic 



limit- the quadratic dependence of the kinetic energy on the ramping time (11) shows 
that the corresponding elastic relaxation occurs on a separated time scale 

1 f2 

For small values of [/, this can be much larger than the first relaxation time. Moreover, 
long ramp-up times delay the onset of elastic relaxation. This is a consequence of 
the reduced deviations of the transient momentum distribution from its equilibrium 
counterpart; hence even less phase space is available for scattering processes and 
thermalization due to Boltzmann dynamics is even more delayed. 

In consequence, the temporal separation of the elastic relaxation by scattering 
events from the earlier energy relaxation spans a transient prethermalization regime. 
Approaching the adiabatic limit, i.e. for long ramp-up times tR or weak interactions 
[/, this time regime becomes more extended in time. This agrees with the physical 
intuition on the momentum relaxation by scattering processes which should be enhanced 
by higher kinetic energies. Such observations parallel those made, for instance, in glasses 
where relaxation times from transient nonequilibrium states grow with decreasing kinetic 
energies (low temperatures). 

The most important observation is that within the prethermalization regime the 
kinetic energy does not equal an effective temperature. Instead, the continued presence 
of an approximate or even sharp [TTj discontinuity at the Fermi energy suggests a zero- 
temperature description for the prethermalized transient state. After thermalization 
has been completed, the kinetic energy determines an effective temperature Te^ defined 



by AE = TT Pi? 7^1/6. Inserting the temperature dependence of AE (11) one 



observes that the effective temperature is proportional to the inverse ramp-up speed. 

Teff ^ — (16) 

7. Consequences for the observability of nonequilibrium BCS dynamics 

As the above results, which have been obtained for the repulsive Hubbard model, depend 
only on even powers of the interaction they can be carried over to the case of an attractive 

§ The linearization of the scattering integral around the zero temperature initial state is valid at the 
onset of thermalization. With increasing temperature, however, the linearization must be performed 
around a Fermi-Dirac distribution of nonzero temperature. This leads to a reduction of the scattering 
integral with time and to the expected asymptotic behavior of the relaxation dynamics. We do not 
follow the Boltzmann dynamics here. 
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interaction [/ < 0. In both cases the initial state is given by the noninter acting Fermi 
gas but the final states after the adiabatic or sudden ramping diflFer. For both attractive 
and repulsive interactions a correlated many-particle state builds up; however, below the 
critical temperature Tc a Fermi liquid with pairwise attractive interactions around the 
Fermi surface exhibits the distinct feature of superconductivity or superfiuidity: Then 
a single particular normal mode of the Fermi liquid becomes unstable with respect to 
Cooper pairing of low energy fermions. 

The pairing behavior can be discussed both within the attractive Hubbard model 
and by referring to the effective BCS Hamiltonian 

HFn?' = mJ2iApi'-^^'^t (17) 

Although both approaches describe pairing in a similar way, they are not fully equivalent: 
the BCS interaction can be considered as that subset of the two-particle Hubbard 
interaction matrix elements which describes the behavior of pairs of fermions with 
opposite momenta. Other two-particle interaction matrix elements of the Hubbard 



model are not represented in the BCS Hamiltonian (17). 

The nonequilibrium behavior of the BCS Hamiltonian (17) following a nonadiabatic 
switch of the interaction A(t) has been recently discussed [20l EH [22l [23l [24]. There it 
has been shown that a quantum quench in the effective pairing interaction A leads to 
the buildup of a paired many-particle state which, furthermore, exhibits characteristic 
nonequilibrium signatures. For example, temporal oscillations of the order parameter 
of the superfiuid transition, the energy gap parameter A(t), have been observed. 

However, experiments in ultracold fermions are well-described by the Hubbard 
model. While it is experimentally feasible to change the relative strength of the Hubbard 
interaction in time, it is not possible to do so only for that subset of matrix elements 
which correspond to the BCS Hamiltonian. Hence the question arises whether the 
dynamics described by the remaining matrix elements interferes with the predicted 
nonequilibrium BCS dynamics. 



7.1. Excitation of nonequilibrium BCS dynamics by finite time ramping of the Hubbard 
interaction 

In the following we will address one aspect of that question: whether the effects of the 
additional heat inserted by the feasible quench of the Hubbard interaction as compared 
to a theoretical quench of the BCS interaction alone already rule out the observation of 
nonequilibrium BCS dynamics. 

Obviously, the inserted energy depends on the ramping time. In order to observe 
nonequilibrium BCS dynamics in ultracold Fermi gases, two competing requirements 
must be met: 

(i) Firstly, the ramping must be sufficiently nonadiabatic to excite nonequilibrium 
behavior. 
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(ii) Secondly, the system must not heat up beyond the critical temperature of 
superconductivity despite the energy intake by the ramping. 

In a previous work [26] it has been pointed out that a ramp-up which is nonadiabatic 
on the scale of the BCS instability, i.e. 

tR < 1/A (18) 

is sufficient to excite nonequilibrium BCS dynamics. As A is typically a small parameter 
and as the kinetic energy depends inversely on the ramp-up time, this eases the problem 
of overheating with respect to a sudden quantum quench. 



7.2. Overheating of BCS system beyond its critical temperature cannot be excluded 
However, plugging the slowest sufficiently nonadiabatic linear ramp tR = tecs =^ 1/A 



(cf. 18) into (16) shows that even for that limiting case the energy intake is not 



negligibly small. Instead, the final temperature is of a similar order of magnitude as 
the characteristic energies involved in the BCS dynamics, namely Ten ^ UA. Hence, 
energetic arguments do not imply a clear separation of energy scales related to the 
nonequilibrium BCS dynamics and temperature, respectively. Therefore, we expect 
that finally the temperature will infiuence and possibly destroy the nonequilibrium BCS 
behavior. 



7.3. Prethermalization opens window for observation of nonequilibrium BCS behavior 

Fortunately, the thermalization of the momentum distribution is delayed by 
prethermalization. For the limiting ramping time tji = 1/A thermalization due to 
elastic scattering processes is deferred to the time scale set by 

t| BCS 1 1 2/pfU ^ , 

For small interaction strength, the elastic relaxation time becomes much larger than 
the time scale of the BCS dynamics. Hence a long-lasting transient time window is 
opened by prethermalization; during this time the momentum distribution around the 
Fermi energy resembles that one of a Fermi liquid at zero temperature and most of the 
excitation energy is virtually stored in high energy modes sufficiently far apart from it. 
Thus we expect the buildup of the BCS instability and the persistence of BCS dynamics 
until quasiparticle scattering redistributes the excitation energies down to low energy 
modes, implying thermalization. 

Based on our observation of prethermalization we are happy to report that 
nonequilibrium BCS dynamics is, in principle, accessible by suitable linear ramping 
protocols performed with ultracold fermions. Other ramping protocols as discussed by 
Eckstein and KoUar can further facilitate these experiments by reducing the effective 
temperature. 
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Figure 6. This sketch compares the relevant relaxation time scales in the case of linear 
finite time ramping on the time scale Ir ~ 1/Abcs (1) with those following a quantum 
quench (2); the later is assumed sudden w.r.t. the time scale set by the hopping 
matrix element. In both cases an initial relaxation (IR) is followed by a transient 
prethermalization regime and final thermalization (Th). The finally achieved effective 
temperature is of the order or larger than the critical temperature of superconductivity. 
Note that the prethermalization regime following a sudden quench may not reach 
into the time regime where signatures of the BCS nonequilibrium dynamics can be 
expected. Hence the long duration of the prethermalization regime in the case of finite 
time ramping is crucial for a clear observability of nonequilibrium BCS dynamics. 

8. Conclusions 

In this paper we have analyzed the implications of linearly ramping up the interaction 
strength of the Hubbard model in more than one dimension and at zero temperature. 
We assumed a finite ramping time and discussed corrections to the adiabatic limit as 
well as the crossover into a quantum quench using second order Keldysh perturbation 
theory in real time. We observed different relaxation behavior for two observables, 
namely the momentum distribution and the kinetic energy. On the one hand, the 
momentum distribution does not relax on a first time scale which is related to second 
order perturbation theory; it shows a characteristic mismatch of the quasiparticle residue 
if compared with adiabatic ramping. This mismatch provides a quantitative measure 
for the nonequilibrium character of the momentum distribution in the first time regime 
after ramping. On the other hand the kinetic energy summed over all momentum 
modes remains unchanged under the subsequent Boltzmann dynamics of the momentum 
distribution; the later becomes effective on a second, much longer time scale and leads 
to the thermalization of the momentum distribution. Consequently, a transient time 
regime characterized by prethermalization appears between these two well-seperated 
time scales. During this regime the excitation energy of the switching process is 
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present predominantly in high energy momentum modes; around the Fermi energy a 
zero-temperature description remains apphcable. For smaU deviations from adiabatic 
switching, i.e. for smaU energy intakes due to the ramping, this regime is particular 
extended in time. We compute the dependence of the final temperature on the ramping 
time and explain how differing expectations have been obtained previously. Finally 
we apply our findings to the question whether nonequilibrium BCS dynamics can be 
observed in ultracold Fermi gases. On the one hand, the final temperature reaches the 
same order of magnitude as the critical temperature of the BCS problem which makes the 
observation of such dynamics in a true long-time limit -at least- very difficult. On the 
other hand, a suitable choice of the ramping speed and the observed prethermalization 
of the momentum distribution function open a long transient time window; during 
that time the excitation energy inserted by the ramping remains located in high energy 
momentum modes. Therefore, nonequilibrium BCS dynamics, which happens on a small 
energy scale around the Fermi surface, is not affected by the inserted energy. We expect 
that prethermalization fosters in a similar way the observability of low-temperature 
nonequilibrium quantum dynamics of other phenomena which also require the presence 
of a sharp Fermi surface. 
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Appendix A. Keldysh perturbation theory for the Hubbard model 

In the following appendix we sketch the real-time Keldysh perturbative approach to 
the Hubbard model. For conventional perturbation theory of the ground state the 
Gell-Mann and Low theorem links the asymptotically free ingoing {t — oc) and 
outgoing {t oc) ground state. Interaction effects are switched on adiabatically, treated 
perturbatively and switched off again on the way of a linear time evolution. Perturbation 
theory beyond the ground state, however, cannot rely on the Gell-Mann and Low 
theorem which only links ground states. Therefore, calculations have to be performed 
on a closed time contour, mapping the outgoing state back into a representation of the 
ingoing states. The Keldysh technique [33j has implemented this contour evolution for 
the calculation of nonequilibrium Greens functions by representing those on a 2 x 2 
Keldysh index space. Hence, Feynman rules have to be set up which include the matrix 
structure of Keldysh Greens functions. Details of the following calculation can be found 
in [31]. 
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Fermionic propagator: 
Boson propagator: 

Boson absorption vertex: 
Boson emission vertex: 



\ Ik'k' 



O (cf. 



Sit - t')S{x - x')bibtl 



/m 

^/2 



V2 



Sk'k" 
^k'k" 



.(1) 



Table Al. Feynman rules for the Hubbard interaction in Keldysh space, n represents 
indices in Keldysh space, a in spin space, r^^)^,, the 2x2 Pauli matrices. The boson 
absorption and emission vertices (which are 3-tensors in Keldysh space) differ in the 
chosen Larkin-Ovchinnikov representation. 



Appendix A.l. Nonequilibrium Feynman rules for the Hubbard interaction 

The tensor structure of Feynman rules for the Hubbard interaction on Keldysh space can 
be constructed in analogy with any exchange interaction: Fundamental vertices describe 
the emission and absorption of exchange bosons by free fermions; working in the Larkin- 
Ovchinnikov representation of Keldysh space, their nontrivial tensor structure can be 



found in [30] and is quoted in table Al A Keldysh representation of the Hubbard 
interaction can be constructed by linking emission (7) and absorption (7) vertices by a 
local, instantaneous boson propagator (D) which is diagonal in spin and Keldysh space. 
The time dependence of the ramping process is included in an explicitly time-dependent 
interaction strength U{t) and symmetrically attributed to both vertices. In the Larkin- 
Ovchinnikov representation the Keldysh matrix Greens function can be composed from 
the retarded (G^), advanced {G^) and Keldysh component {G^) Greens function. 

/ qR qK \ 

G^^>{x,t,a^,x\t',G'^)= \ ^ \{x,t,G^,x',t',G'^) (A.l) 

For noninteracting fermions, G^, G^ and G^ are given by 

g5(x, t, a., x\ t', a;) = -^Y: e^'^(— )--(*-*^) ( ^^'^ " T ^^'^ ] (A.2) 



This makes the free fermion propagator explicit. Diagrams can be evaluated using real- 
space Feynman rules for Greens functions. Note that whenever propagators are linked 
to vertices, the related Keldysh indices have to be contracted. 



Appendix A.2. Second order Keldysh perturbation theory for the Hubbard model 



Fortunately, perturbation theory for the Hubbard model is simplified by various 
symmetries: the bare first-order Hartree diagram is cancelled by a symmetric definition 



Crossover from adiabatic to sudden interaction quenches in the Hubbard model 22 




Figure Al. The only diagram contributing to a second order perturbative expansion 
of the Keldysh Greens function. 



of the Hubbard interaction, renormalized Hartree diagrams vanish to aU orders because 
of particle-hole symmetry and Fock diagrams do not arise because the interaction is 
effective only between particles of different spin. In consequence, only a single diagram 
contributes to a second order perturbative calculation of the Keldysh Greens function 



(cf. figAl). 



Appendix A. 2.1. Polarization operator. In order to make the explicit expressions more 
handy it is advisable to split off the calculation of the closed fermionic loop, known as 
the polarization operator 11. Moreover, the shorthand notation i = {xi^ti^a^} where i 
is an integer number is introduced. To avoid confusion, it does not extend to Keldysh 
indices. The later are, generically, denoted by Greek letters or, explicitly, by their values 
(1 or 2). 



n.3..(3,2)= %%3(*3)G,3.4(3,2)G.3.,(2,3)7,%,(t2) (A.3) 

Performing the contractions of the Keldysh indices allows to represent the polarization 
operator as a 2 x 2 matrix in Keldysh space. 



n(3,2) = 



I G^(3,2)G^(2,3) 
+G^(3,2)G-^(2,3) 

V 



G^(3, 2)^^(2, 3) + 6*^(3, 2)^^(2, 3) \ 
+G^(3,2)G^(2,3) 



G^(3,2)G^(2,3) + G^(3,2)G^(2,3) / 

(A.4) 



This defines naturally the components 11^ = IIn, 11^ = 11x2 and 11^ = 1122 



Appendix A. 2. 2. The ^setting sun^ diagram. The second order correction to the Keldysh 
Greens function follows from the evaluation of the diagram figure Al which can be 
written as 

G21^,(1,10 = ^ / dx2dx^dx^dx^ J dt2dt^dt^dt^ ^ ^t.^mO^^'^hm^^i^^) 

X D,3,, (3, 5)n.,«, (5, 4)D«,,, (4, 2)G,,,, (3, 2)r' (^2)6,,,, (2, 1') (A.5) 



To achieve a second order correction G denotes on the right hand side the 
respective noninteracting Greens function (superscripts are suppressed) . Performing all 
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contractions over Keldysh indices and the integrations over trivial bosonic propagators 
leads to 





(l,3)G«(3,2)n«(3,2)G2M(2,l')+ 


(A.6a) 




(l,3)G^(3,2)n^(3,2)Gi,,(2,l')+ 


(A.6b) 




(l,3)G^(3,2)n^(3,2)Gi,,(2,l')+ 


(A.6c) 


G^il 


(l,3)G^(3,2)n^(3,2)G2.,(2,l')+ 


(A.6d) 


G^i2 


(l,3)G^(3,2)n^(3,2)Gi,,(2,l')+ 


(A.6e) 


G^i2 


(l,3)G^(3,2)n^(3,2)G2.,(2,l')+ 


(A.6f) 




(l,3)G«(3,2)n^(3,2)Gi,,(2,l')+ 


(A.6g) 


GkiI 


(l,3)G^(3,2)n^(3,2)G2.,(2,l')+ 


(A.6h) 


G^i2 


(l,3)G^(3,2)n^(3,2)G2«,(2,l')} 


(A.6i) 



From this result second order corrections to all components of the Keldysh matrix Greens 
function can be extracted by specifying the remaining Keldysh indices /^i and 



Appendix A. 3. Second order correction for the Keldysh component Greens function 

Evaluating the Keldysh component Greens function is sufficient to extract the 
information on the momentum distribution function (cf. [5]). This is done by setting 
/^i = 1 and = 2 in (A. 6). Inserting (A. 2) together with (A. 4) into the Keldysh 
component of (A. 6) shows that - because of mismatching time restrictions imposed by 



the involved B-functions - the terms (A. 6c) and ( A.6e) vanish. The remaining terms can 



be paired such that all pairs share the same momentum sums and spatial dependences. 



Gg),^,(l,l') = dx^dx^ j dhdt, 



^/U{h)U{h) 



G^(l, 3)G^(3, 2)n^(3, 2)G^(2, l') + G^(l, 3)G^(3, 2)n^(3, 2)G^(2, l')+ ([A^l+S) 
G^(l, 3)G^(3, 2)n^(3, 2)G^(2, 1') + G^(l, 3)G^(3, 2)n"^(3, 2)G^(2, 1')+ (1X712+4) 
G^(l, 3)G^(3, 2)n^(3, 2)G^(2, 1') + G^(l, 3)G^(3, 2)n^(3, 2)G^(2, 1') (|Ajl5+7) 



G^(l, 3)G^(3, 2)n^(3, 2)G^(2, 1') 



dAjje) 



Although they also share a common time dependent phase factor e^(^2-^3)(ep-ep/+eg2-egi)^ 
their time dependencies %i and phase space factors vary. Integrating out the internal 
positions makes momentum conservation explicit. Then the Keldysh component Greens 
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function can be written as a sum over four pairs 



G^(2)(l,l') = >^ E 



'X 



p'pqiq2 
i=l..A 



%i{ti,tr,Ae)^i{npr,np, Uq^ , nq^)5{p' - p - q2 + qi) (A. 



The phase space factors can be easily read off 

^i+3{np',rip,nq^,nq^) =- {np> - Up) 
^2-rA{np'^np,nq^,nq^) = (1 - 2n^J(l - 2nq^){np' - Up) 
^5+7(v^ n^i, n^J = (1 - 2n^J[np/ + - 2np>np] 
^^{np'.np, nq^.Uq^) = - {1 - 2nq^)[np> + rip - 2npmp 



(A.9,1+3) 



(A.9,2+4) 



(|A.9|5+7) 

The time dependent integration kernels are evaluated for the linear ramp-up scenario 
Q, assuming that ti^ty > tR. This means, in a strict interpretation, that the obtained 
result is only applicable to study the behavior of the momentum distribution function as 
soon as the full interaction has been reached. Transient phenomena during the ramp-up 
are neglected. Then integrating out internal times leads to 

1 2(1 -cos(A6ti^)) 



1+3 



2+4 



5+7 



f2 



+ 



1 



T (Ae)3 
1 



(A.IO) 

tv — ti 



+ 1 



[(Ae)2 ■ " Ae 
1 2(1 - cos(AetH)) 



f2 



{Aey 



+ 



1 



tR (Ae)3 
1 



JAe(tR-ti) 



(A.ll) 



(Ae)^ 



iAe(ti/-ti) 



Obviously, this calculation has produced a secular term in (A.IO) which is proportional 
to the difference of external times ty — ti. However, this term does not affect the 
momentum distribution. 



Appendix A. 4- Second order correction to the momentum distribution 

Since the momentum distribution follows from the Keldysh component Greens function 
at equal times ^ , the secular terms in its second order correction drop out from the final 
result. Moreover, this implies that for linear ramping procedures and the assumption 
ty = ti all time dependent integration kernels T1+3 = = 2^5+7 = =: T coincide. 

Factorizing 3 from the second sum in (A. 8) allows to add up the phase space factors 
of the paired terms. In total, this generates 

^ = -4[n^^np/(l - np){l - Uq^) - (1 - n^J(l - np')npnq^] (A.12) 
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Inserting J dE S{ep — epf + eq^ —E) and canceling prefactors allows to define the fermionic 
phase space factor 



p'pq2 



(A.13) 

Then reading oS the Fourier transform we arrive at the second order correction to the 
momentum distribution given in 



Appendix B. Review of properties of the phase space factor Jk{E,n) 



The fermionic phase space factor Jk{E^n) has been discussed previously [34l|35]. Here 
we review some of its properties. Firstly, it shows quadratic asymptotical behavior for 
small energies around the Fermi energy ep- 



(B.l) 



Then we denote a second proportionality with B(x) being the Heaviside step function 
and Ik{E) a suitable reduced phase space factor. 



Jk{E, n) - [Qi-eumE] - e(efc)e(-E)] hm) ^ 



Jk{E <0) = Vefc<0 
JkiE > 0) = Vefc > 

(B.2) 

This implies that the phase space factor cuts off all formal energy divergences in ([g]). 
Thirdly, we evaluate the reduced phase space factor in the limit of infinite dimensions. 



In this limit all three momentum sums in (A.13) are promoted to energy integrations 



with respect to a Gaussian density of states ([2]). This already includes momentum 
conservation [28] . While the first integration is trivial (energy conservation) , the second 
one can be performed analytically. This gives 



h{E) 



V2wh* 



, * ^ -—■ 
h — e 4h* 

2 



/' 

Jo 



de e 4h*2 



erf 



h* 



erf 



E 



h* 



(B.3) 



A first order expansion of the error function effectively decouples two internal energy 
integrations such that an analytic approximation becomes possible 



Ik{E) 



1 



E^ e 6h*2 



(B.4) 



Fig Bl compares this approximation with a full numerical evaluation of (B.3). As for 



all numerical computations, the bandwidth has been set to six times the hopping rate 

and an energy resolution of A = 6000 sites has been used. 
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Figure Bl. Comparison of the analytical approximate calculation (B.4) and a full 
numerical evaluation of the phase space factor Ik{E). Prefactors are suppressed. One 
observes that the approximation overestimates the available phase space. 
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